High Spatial Resolution Spectroscopy of W51 IRS2E and IRS2W: 
Two Very Massive Young Stars in Early Formation Stages 
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ABSTRACT 

We present i^-band spectra of the near infrared counterparts to IRS2E and 
IRS2W which is associated with the ultracompact HII region W51d, both of 
them embedded sources in the Galactic compact HII region W51 IRS2. The high 
spatial resolution observations were obtained with the laser guide star facility and 
Near infrared Integral Field Spectrograph (NIFS) mounted at the Gemini North 
observatory. The spectrum of the ionizing source of W51d shows the photospheric 
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features NIII (21155 A) in emission and Hell (21897 A) in absorption which lead 
us to classify it as an young 03 type star. We detected CO overtone in emission 
at 23000 A in the spectrum of IRS2E, suggesting that it is a massive young object 
still surrounded by an accretion disc, probably transitioning from the hot core 
phase to an ultracompact HII region. 

Subject headings: stars: early type — stars: formation — infrared: stars - 
techniques: spectroscopy 



Introduction 



The mechanisms of massive star formation (M* > 1OM ) are the subject of intense 
discussion even after two decades of great progress and improvement of observational tech- 
niques. In the last two decades, the development of larger and more sensitive infrared 
detectors associated with the commissioning of 8-m class telescopes equipped with adaptive 
optics has made it possible to observe ever deeper into the sites of massive star formation at 
relatively short wavelengths. Yet, while all these improvements have contributed to our un- 
derstanding of massive star formation, the basic process, i.e., formation via an accretion disk, 
coalescence, competit ive accretion, or a comb ination of these is still unknown (for a review 
on this discussion see IZinnecker fc Yorkd 120071 ). The Spitzer Space Telescope is providing an 
unprecedented view at longer wavelengths and soon the Herschel Space Observatory and 
ALMA will combine very high sensitivity and angular resolution to probe the earliest phases 
of massive star formation. 

W5 1 is one of the most luminous co mplexes of massive star for ming regions in the 
Galaxy (IGoldader fc Wynn- Williams! Il994l ) with multiple HII regions (j Wilson et al.l Il970l ) 
that host at least s ix HII regions with embedded clusters, all of them optically obscured 
(IKumar et al.l I2004J ) . The most luminous of all HII regions is G49.5-0.4 which can be re- 
solved into mainly two bright infrared sourc es: IRS1 and IRS2, which are coincident with 
radio sources W51e and W51d respectively (jMartinl Il972l ). IRS2 is a compact HII region 



harb oring three ultracompact HII (UCH II) regions W51d, W51di and W51d_2 (iMehringer 



2001 



1994) and a plethora of infrared sources ( 



Okamoto et al. 2001: Kramer et al 



Goldader fc Wvnn-Wi 



2001 



Figueredo et al 



lams 



2008 



1994 



Okumura et al 



Barbosa et al.l 120081 ). 



The most prominent sources in the K band are named IRS2 East and IRS2 West (or IRS2E 
and IRS2W, respectively, for short). IRS2W wa s ide ntified as a NIR a nd M IR peak of nebu- 
lar emission (IGoldader fc Wynn- Williams! Il994l ) and lOkamoto et al.l (120011 ) have correlated 
it with the UCHII region W51d. 
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In this letter we report the results of spectroscopic observations taken with laser guide 
star (LGS) adaptive optics of IRS2E and the UCHII region W51d (IRS2W). In the section M 
we report the technical details of the observations and in the section [3] we report and discuss 
the results. Finally, in the section [5] we summarize our conclusions. 

2. Observations 

The data were obtained with the Near infrared Integral Field Spectrograph (NIFS) 
mounted on the Frederick C. Gillett 8-m telescope at the Gemini North observatory on 
Mauna Kea, HI, USA, in queue mode on May, 4th 2007 under program ID GN-2007A-Q-34. 
NIFS provides 2D integral field spectroscopy with R~5200 over a 3"x3" field of view with 
0.043x0.10" rectangular "pixels". Also, we used the Gemini North Adaptive Optics system 
Altair with a LGS to achieve an angular resolution of <0.2" measured upon the 2D images 
build from the datacubes. 

The observations were carried out with the standard procedures for the near infrared 
(NIR): we observed IRS2 in a series of three pointings of 1000 seconds on source and 1000 
seconds on a "blank" sky position at ~ 10" north of IRS2 for each field. The latter image 
was used to subtract the sky emission from the images of the targets. We also observed the 
A0V star HIP 98640 with similar airmass as ISRS2 in order to remove the telluric absorption 
lines from the spectra of the observed targets. The spectra of A0V stars show only the Br7 
feature in absorption in the K band. We eliminated this feature by fitting a Voight profile 
to the Br7 leaving a spectrum with only telluric absorption features. This telluric spectrum 
was used to correct the spectra of the targets. 

Wavelength calibration was achieved with an exposure of an Argon arc lamp obtained 
with the Gemini calibration module GCAL. All data were precessed under IRAe3 envi- 
ronment with scripts written for NIFS available through the Gemini IRAF package. We 
extracted the spectra through a 0.4" circular aperture centered on the sources. This aper- 
ture was chosen to minimize the contribution of the nebular emission from the HII region. 

We show in Figure [1] the finding chart of IRS2. It is composed by the high resolution 
image taken with the adaptive optics NIR camera NACO in the K band. This image is 
available at the ESO public database. We overplotted the radio continuum emission at 2 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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cm from IWood fc Churchwelll (119891 ) in white contours. The black contours represent the 
MIR emission at 12 /xm from the Gemini South MIR camera T-ReCS. All NIFS pointings 
are represented by white squares and the sources which we have spectra are labeled. In 
the following section we present the results of the observations of IRS2E and the UCHII 
region W51d in the "pseudo long-slit" mode extracted from the images delivered by NIFS. A 
complete analysis of these data, including the spectra of all observed sources and the results 
from the IFU mode, such as the velocity maps of the regions observed will be presented in 
a forthcoming paper (IBarbosa et al.l 120081 ). 



3. Results and Discussion 



3.1. IRS2E 



Figure [2] shows the spectrum of IRS2E, the brightest K— band source in the compact 
HII region. It was obtained from the first pointing of the telescope. The continuum shows a 
steep, rising slope to longer wavelengths reflecting the fact that this source is under a heavy 
extinction: Ay ~63 mag according to our unpublished MIR data. This high extinction 
makes IRS2E undetected at wavelengths shorter than 16000 AB The spectrum of Figure 
[2] shows strong lines of Hel (20587 A) and unresolved Br7 (21661 A) that may arise from 
the circumstellar gas heated by hot stars. We also note nebular lines of Hel (21127 A) and 
the H 2 (1-0 (SI), 211 28 A) that come from the diffuse nebular em ission of the HII region. 



Previous authors (e.g. lOkumura et al.ll200ll ; IFigueredo et al.l 120081 ) describe the presence of 



[Felll] lines in the s pectrum of IRS2E. These emission l ines are typical of shocked gas in the 
vicinity of MYSOs (IHanson et al.ll2002l ; iBik et al.ll2006l ). The spectrum in the Figure [2] does 



not show any of these lines. This is due to the small aperture used to extract the spectrum. 
Actually, our unpublished line maps show that [Felll] emission comes from the extended 
region between IRS2E and OKYM2/KJD5. 

More interesting, IRS2E shows CO overtone emission as we detected the (2-0), (3-1), 
(4-2) and ( 5-3) bandheads between 22935 A and 2383 8 A. Previous NIR spectroscopy of 



IRS2E (e.g.. lOkumura et al.ll200ll ; IFigueredo et al.ll2008l ) could not detect the CO bandhead 



emission due to the large contribution of nebular emission through wide and long slits. 

The detection of the CO overton e in emission is a ssociated with the pres ence of an 
accretion disc, both in low-mass (e.g., iNajita et al.l 119961 ) and high-mass YSOs (IBik fc Thi 



2 Unpublished images obtained from the AO assisted camera NACO at the VLT shows that the source 
detected at the J band by previous authors is just a bright knot of gas. 
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2004 ; iBlum et al.ll2004j ). The CO overtone emission comes from the inner parts a circum- 



stellar disk, in regions of high column densities (10 20 -10 21 cm 2 ) and temperatures between 
1500-4500 K. 

IRS2E is a deeply embedded source which is bright at wavelengths longer than 20000 A 
and is not associated with any UCHII region. These characteristics are those expected for 
a very young and massive YSO that is transitioning to an UCHII region; the high accretion 
rate presumably prevents the formation of an UCHII region as the ionized gas near th e sourc e 
falls onto the star, according to the evolutionary scenario proposed by IChurchwelll (120001 ) . 
This class of object was observed in NGC 3576 by iBarbosa et al.l (120031 ). for example. 



3.2. W51d 



Figure E] shows the spectrum of the UCHII re g ion W51d. Treat e d as a peak of nebulosity 
in the NIR bv lGoldader fc Wvnn- Williams! (119941 ). iFigueredo et all (120081 ) identified a stellar 
object as responsible for the ionization of this region based o n the correlation between the 
radio and MIR emission presented by lOkamoto et al.l (120011 ) and the new high resolution 
images obtained by us or collected at public databases. 

The spectrum of W51d was normalized and zoomed near the intense Br7 line to better 
show the detected features. This line has a broad component (~ 700 km s^ 1 ), probably 
originated in fast winds coming from the hot ionizing star and two scenarios are possible 
to explain it. In the first case, the very young star has strong winds that blew away its 
surrounding dust. In the second case, the ionizing star is still surrounded by a torus that is 
seen more face on. However, only higher resolution data could favor one of these scenarios. 
Besides this nebular line, we detected photospheric lines corresponding to the NIII multiplet 
near 21155 A in emission and Hell at 21897 A in absorption. 

The ionizing source of W51d can be classified as an early O star, given the detected 
features in the spectrum of Figure Co mparing this spectru m with the catalog of spectra 
of optically visible OB stars presented by lHanson et al.l (120051 ) we can constrain its spectrum 
type to an 03-4 type star. The spectrum of an optically visible 03V star shows both the 
NIII multiplet (in emission) and the Hell (in absorption) at wavelengths longer than 20000 
A. Although our spectrum cannot establish its spectral subtype, the absence of CIV (20780 
A) in emission, as seen in mid type O stars, favors the identification of this object as a star 
hotter and more massive than type ~ 04. It would make this object one of the most massive 



ioniz i ng source of an UCHII re gion (I Watson fc Hansonl Il997l ; lHanson et al.l |2002| ; iBik et al. 
20051 : IBlum fc McGregor! l2008h . 
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Wood fc Churchwelll (119891 ) and iMehringerl (119941 ) classify the ionizing source of W51d 



as an 05.5 and 05 respectively, both based on radio continuum emission assuming a distance 
of 7 kpc to this source. The disagreement between the results based on radio data and the 
spectroscopy can be explained if the massive YSO had enough time to clear the circumstellar 
material, lowering its density (and the extinction as well). In this case, the star will not be 
as embedded as a hot core (like IRS2E) and some fraction of the total number of the Lyman 
continuum photons emitted by the star will freely escape without ionizing the surrounding 
gas. If this is the case, the ionizing source would be of earlier spectral type. 



3.3. The Distance to IRS2 

Normally, one could derive the distance to W51d using its photometry in the NIR and 
the spectral type obtained above. However, we could not obtain a reliable photometry in 
the H band from the NACO image. W51d seems to be saturated (as well as the object 
BBCDF4) and every attempt to reconstruct its PSF has failed. Moreover, this source seems 
to have strong excess emission and in this case we cannot derive its H — K color nor assign 
its extinction. Without its H — K color and extinction we cannot derive an estimate for its 
distance using spectroscopic parallax. 

On the other hand, based on the spectral type presented in the previous section, we 
can estimate the distance to IRS2 using the radio data available in the literature. We can 
derive the number of photons emitted in the Lyman continuum by an 03 star and scale the 
distance until we reproduce the flux measured in the radio at a given wavelength. We need 
to keep in mind the limitations of this method, since it assumes that all photons emitted in 
the Lyman continuum are effectively ionizing the gas, and this may not be true as explained 
in the previous section. Some amount of the total number of photons may be escaping the 
UCHII while other photons are destroyed by dust. 



Starting with equation (4) given by iKurtz et al.l (119941 ) (£ = 1, i.e. all photons in the 



Lyman continuum effectively ionize the gas) we obtain the distance d (in kpc) as a factor 
of the electron temperature T e , the correction factor a(is,T e ), the frequency v and the flux 
density S v as: 



The critical parameter in the above equation is the number of photons emitted in the 
Lyman continuum Nc- Values for that parameter can vary an order of magnitude depending 
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on which author is quoted. We derived Nc as follow: for spectral types ranging from an 03 
to an 04 star at the zero age main sequence (ZAMS) we took their effective temperature 
from Table 1 of lMartins etaD (l2005h: 4 4.616 K for an 03 and 43,419 K for an 04. Using the 
Geneva models of ISchaller et al.l (119921 ) we obtained the bolometric magnitude M^i for each 
star in the ZAMS (-8.94, 03 and -8.7, 04). We derived the stellar radius from the effective 
temperature above and its corresponding luminosity: 5.476 and 9.2 R & for an 03 and 4.638 
and 8.7 R Q for an 04. The radius and the luminosity calculated for a star in the ZAMS are 
significantly smaller than those quoted for main sequence stars and may be more reliable for 
W51d, given its young age. Finally, we used the radius for each star to obtain the number 
of ionizing photons emitted per unit time from the equation (9) and the ionizing fluxes Q n 
for each respective effective temperature also quoted in the Table 1 of lMartins et al.l (120051 ). 



According to these procedures an 03 can produce 1.86 x 10 49 s 1 and an 04, 1.44 x 10 49 s 1 
in the ZAMS. 

Returni ng to our equation (1), S„ and T e quoted for W5 1d are: 5.35 Jy and 10 4 K a t 
v =15 GHz JWood fc Churchwelj[l989T ) and a{u,T e ) = 0.9767 iMezger fc Henderson! Jl967h 
The r adio flux of W51d is blended with that of W51di in data presented by iMehringer 
( 11994J ). so it is overestimated and therefore cannot be used. The results are: for a ZAMS 
03 star d = 5.8 kpc and for an ZAMS 04 d — 5.1 kpc. We can compare o ur res ults with 
those obtained using the number of ionizing photons given by lMartins et al.l (120051 ) for main 
sequence stars, in this case the distance is 8.8 kpc for an 03V and 7.4 kpc for an 04V. 

We cannot take these results at their face value, given the limitations of the method 
described above. However, the results can be used as upper limits to the distance to IRS2. 
The distance to W51 North, which may be associated with IRS2, quoted in the literature 
were obtaine d through observatio ns of proper motions of water masers and amounts to 8.3 
(± 2.5) kpc (jSchneps et al.lll98ll ) and 6.1 (± 1.3) kpc (llmai et al.l 120021 ) . The upper limits 
obtained above are compatible (within the errors) with those obtained from the water masers 
if we consider the ionizing star of W51d an 03-04 star in the ZAMS. However, these upper 
limits would make IRS2 part of W51 Main just as an effect o f projection on the sky , given 
the distance of 2 kpc obtained with spectroscopic parallax by iFigueredo et al.l (120081 ). 



4. Conclusions 

We presented K-band spectra of the embedded sources IRS2E and IRS2W/W51d in the 
compact HII region IRS2. The high angular resolution achieved by NIFS and Altair made 
possible the separation of the strong extended emission from the emission of the MYSO 
itself. It has proved a valuable tool in the study of MYSOs as we summarize our results 
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presented in the previous sections as follow: 

1) We identified the ionizing source of the UCH1I region W51d as a massive star ear- 
lier than about type 04 based on the NIII multiplet (in emission), the Hell line (in 
absorption) and the absence of CIV in emission; 

2) We detected compact (unresolved) CO overtone emission toward source IRS2E for the 
first time, suggesting the presence on an accretion disc around this source. The absence 
of radio emission and the presence of an accretion disk indicate that IRS2E is a massive 
YSO. 

3) Combining the spectroscopy of W51d and its radio data, we found an upper limit to 
the distance to IRS2 between 5.1-5.8 kpc for an ionizing star of spectral type between 
03-04 in the ZAMS. 

The fact that we detected two very massive stars, close in space, one which has cleared 
its surrounding material enough to probe its photosphere and another in an younger phase, 
deeply embedded in its dust cocoon and still surrounded by a disk, indicates that the evolu- 
tive timescale of MYSOs is extremely short. 
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Fig. 1. — Finding chart of IRS2. The white contours are the radio continuum emission at 2 
cm. The black contours represent the MIR emission at 12 ^m. White squares represent the 
field of view of NIFS (3"x3") and indicate the positions of the three pointings of our program. 
The sources labeled are those which we have "pseudo long-slit" spectra. All contours are in 
arbitrary units. 
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Fig. 2. — K-band spectrum of source IRS2E. After the spectrum was extracted it was 
multiplied by a black body of 9000 K to recover the true spectral shape of the source. The 
flux is in arbitrary units. 




Fig. 3. — fT-band spectrum of source W51d. This spectrum was zoomed to show the NIII 
multiplet and Hell. The thin line near 20587 A is a residual feature produced after the 
correction of the telluric lines. 



